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Abstract 



o 



We analyse photoproduction of r/ mesons off the proton in a gauge-invariant chiral unitary framework. The interaction kernel for 
meson-baryon scattering is derived from the leading order chiral effective Lagrangian and iterated in a Bethe-Salpeter equation. 
The recent precise threshold data from the Crystal Ball at MAMI can be described rather well and the complex pole corresponding 
to the S n(1535) is extracted. An extension of the kernel is also discussed. 
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Introduction: Coupled-channel unitary extensions of chiral 
perturbation theory have been established as a viable tool to 
investigate the chiral SU(3) dynamics of QCD, for early pa- 
pers see e.g. IE B B S] • The main reason for this is that 
within such a framework, resonances close to and even below 
the relevant thresholds can be generated dynamically. The two 
premier examples are the A(1405) that features prominently in 
antikaon-proton scattering and the S n (1535) that dominates the 
1 threshold cross section in eta photoproduction off protons. In 
the context of unitarized chiral perturbation theory, this was 
first investigated in Refs. @, 0] based on the next- to-leading 
order (NLO) effective Lagrangian (i.e. using the pertinent con- 
tact interactions). The inclusion of p- waves was studied in 
(SI], and the role of the titiN final-state was investigated in 
13]. T] and K]' photo- and electroproduction based on a U(3) 
extension of the chiral Lagrangian was considered in 11011 . In 
these ground-breaking papers, the issue of gauge invariance was 
jnot considered. In Ref. [11], a gauge-invariant framework for 
meson photoproduction was developed and applied to associ- 
ated strangeness production, using the leading order Weinberg- 
Tomozawa interaction as the driving term in the Bethe-Salpeter 
equation. In this note, we extend this method to the process 
yp — > rjp and in particular to a determination of the mass and 
width of the 5 n (1535). Eventually, NLO contributions will 
have to be included, a first step for the purely hadronic interac- 
tions within our approach was recently reported in 11211 . 

Framework: In this work we follow the construction proce- 
dure of the minimal approach to meson photoproduction, which 
is exactly unitary in the subspace of meson-baryon channels 
and gauge-invariant, as developed in Ref. [11]. We refer the 
reader to that paper for a more detailed discussion. The first 
building block of the photoproduction amplitude is the meson- 
baryon interaction, for which we consider the chiral effective 
Lagrangian of QCD at leading order: 



where (. . .) denotes the trace in flavor space, D^B = d^B + 
[T^B], mo is the baryon octet mass in the chiral SU(3) limit, 
while D and F are the axial coupling constants. The relevant 
degrees of freedom are the Goldstone bosons described by the 
traceless meson matrix U = u 2 = exp(/0/Fo), where Fq is the 
meson decay constant in the chiral limit. The meson- and the 
low-lying baryon-fields are collected in traceless matrices cp and 
B, respectively. Moreover, we use = iu\d^U - i[v^, U])u^ 
the so-called chiral vielbein. The quark charge matrix is Q = 
diag(2/3, -1/3,-1/3), and the external vector field is encoded 
in = -eQN 1 . 

The expansion of the chiral connection = [u^d^u]/! - 
iiu^v^u + wv /i w t )/2 in meson fields leads to the meson-baryon 
vertex of the leading chiral order, the Weinberg-Tomozawa 
(WT) term. However, at first chiral order, there are also the 
Born graphs, describing the ^-channel and w-channel exchanges 
of an intermediate nucleon. The full inclusion of these graphs 
in the driving term of the Bethe-Salpeter equation (BSE) leads 
to conceptional and practical difficulties, which have not yet 



4* = ~ m o)#> + -^-{By,y 5 [u», B] ± ) , (1) 
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been solved to the best of our knowledge, see 11121 1 1 3L 11411 
for a more detailed discussion of this issue. However, most 
chiral unitary approaches restrict their meson-baryon potential 
to this interaction, which generates the leading contribution to 
the s-wave scattering lengths. This approach has been remark- 
ably successful in many cases, thus we iterate the WT-potential 
Vwt = g(ii + <f 2 ) v i a me BSE in d space-time dimensions to 
infinite order as follows 

T(f 2 4 1 ;p) = V WT (t 29 t 1 )+ 

J ^VwT^JWtp-l^WiUq^pX (2) 

where the in-, outgoing meson and the overall four-momentum 
are denoted by q\, q 2 and p, respectively. The baryon 
and the meson propagators are represented by iS(fi) = 
-m + ie) and iA(k) = i/(k 2 - M 2 + ie), in order. More- 
over, every element of the above equation is a matrix in the 
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channel space. For the case of 77-photoproduction consid- 
ered here, we restrict ourselves to the following channels: 
{pn , nn + , prj, AK + , lPK + , £ + ^ }. 

Because the kernel of the integral equation Eq. © stems 
from the contact interactions only, we have to deal at most with 
one-meson-one-baryon loops, which are all dimensionally reg- 
ularized throughout this work. It is not possible to express the 
terms necessary to absorb the divergencies in the BSE as coun- 
terterms derived from a local Lagrangian. However, it is pos- 
sible to alter the loop integrals in the solution of the BSE in a 
way that is in principle equivalent to a proper modification of 
the chiral potential itself, see H 1 311 . In this spirit we apply the 
usual MS subtraction scheme, keeping in mind that the modi- 
fied loop integrals are still scale-dependent. This regularization 
scale (p) is used as a fitting parameter, reflecting the influence 
of higher order terms not included in our potential. 

The functional form of the driving term allows to construct 
an explicit solution of the BSE, see fllll . Starting from the cor- 
responding scattering amplitude T((f 2 , p), where exact two- 
body unitarity is guaranteed by construction in the BSE frame- 
work, and following the recipe of fllll we are now able to con- 
struct the gauge-invariant photoproduction amplitude in a most 
natural way, without any use of "artificial restoration" of gaug e 
invariance, i.e. adding contact interactions, see e.g. 
In our approach, we simply couple the photon to any in- and 
external line as well as to the (momentum-dependent) vertices. 

As an intermediate step we construct the amplitude T for the 
process p — > Bcp starting from the potential derived from the 
leading order Lagrangian Eq. Q) V = fYsg, where q is the 
outgoing meson momentum and g is the (D-, F-dependent) 
coupling constant corresponding to the second term of Eq. (Q]). 
Now we add the loop contribution that accounts for the final 
state interaction as follows 

r d d l 

r(j>,0 = v(f> + J — d T^j\ P )is{p-immi). 

Note that the process in question is meson photoproduction off 
the proton, thus T, V and g are 6-vectors in the channel space. 

In the next step we couple the photon in every possible place 
to the hadronic skeleton. A photon can couple via: Bcpy —> (pB, 
arising from the chiral connection, By — > B, (py — > or via 
the Kroll-Rudermann interaction yB — » B(p, stemming from 
the chiral vielbein. These interactions give rise to nine dif- 
ferent topologies presented in Fig. [T] Following the conven- 
tions of 111711 the most general Lorentz covariant form of the 
amplitude for the process: y(k)Bi(p - k) — > Bf(p - q)(pf(q) 
can be written as Tp = ie^Uf{^ i=x SiN^)ut, where is the 
photon polarization vector and Ujj are the initital and final 
dirac spinors, respectively, which are normalized like uu = 
2m. Moreover & denote the coefficients of the pseudo-vectors 
N" e {y 5 yt, 2y 5 P», 2y 5 q», 2y 5 k», y 5 y, y 5 ^, 75**?, 7sW 
with P = 2p - q - k. 

For further study we fix the axis of quantization to the z- 
axis and rewrite the amplitude once more in terms of CGLN 
amplitudes and Pauli spinors and martices as 

T fi = %7t^~sx\Txu (3) 



Figure 1 : Different topologies contributing to gauge-invariant photomeson pro- 
duction off the proton. Solid, dashed and wiggly lines denote baryons, mesons 
and photons, in order. The filled circles/squares denote the meson produc- 
tion/the meson-baryon interaction, i.e. T and T, respectively. 

where T = i(cr • e)T\ + (cr • q)(cr • [k X e]) Ti + i(cr • k)(q • e)T^ + 
i(cr -q)(q- e)T^ + i(cr -k)(k- e)T$ + i(cr -q)(h e)T(> + i(cr -q)(k- e)T^ - 
i(cr • q)€{)Ti - i(cr • k)eoTg with k and q normalized three- vectors. 
For the exact form of CGLN amplitudes in form of coefficients 
& we refer the reader to the 111 111 . Let us note here that due to 
current conservation two of the eight CGLN amplitudes can be 
expressed in terms of other six. Moreover two of the remaining 
six CGLN amplitudes are accompanied by scalar components 
of 6 only and thus have no influence on process including real 
photons. In view on photoproduction this leaves us with four 
independent CGLN amplitudes. 

Let us also mention that by construction of the unitary 
hadronic interaction the photoproduction amplitude obeys the 
requirement of two-body unitarity in the subspace of meson- 
baryon channels automatically. There are five different unitar- 
ity classes, which obey the two-body unitarity by themselves. 
Gauge invariance is, however, fulfilled only if all topologies are 
taken into account. Note that within the approximation used 
here, crossing symmetry is violated. 

Results: We are now able to confront our model with the ex- 
perimental results, for which we consider the recent measure- 
ment of the differential cross sections from the Crystal Ball at 
MAMI 0]. The parameters of the model are the unknown 
renormalization scales, whose number is restricted to four due 
to isospin-symmetry: {yu^, l^rjN, J^ka, Hkt}- In contrast to fllll 
we do not consider the meson decay constants as free param- 
eters, but fix them to their physical values. Our input parame- 
ters are: F n = = 0.0924, F K = 0.113, M n o = 0.135, 
M n+ = 0.1396, M v = 0.5478, M K+ = 0.4937, M K o = 0.4977, 
m p = 0.9383, m n = 0.9396, m A = 1.1157, m s o = 1.1926, 
m s+ = 1.1894, D = 0.8 and F = 0.46 (all masses and decay 
constants in units of GeV). 

The p- and d-waves become important already at moder- 
ate energies away from the threshold. Since the kernel of the 
BSE considered here is the WT-term only, which mainly pro- 
duces an s-wave contribution, we restrict the center-of-mass 
(cms) energy of the present analysis to 1487.8 MeV < E cms < 
1609.9 MeV. Note that the WT interaction kernel produces 
a small p-wave contribution via the lower components of the 
Dirac spinors, however the explicit angular dependence of the 
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Figure 2: Best fit of our model (solid lines with shaded error bands) compared to the experimental data from the Crystal Ball at MAMI (filled circles with error 
bars). The differential cross section dcr/dQ, in \jibj sr] is plotted versus cos(#) for the cms energies given on top of each panel. 



amplitude is only described when higher order potentials are 
taken into account |8]. The calculation of the BSE with the 
potential at next- to-leading order is performed in Ref. 1 12] for 
the hadronic part and is to be extended to photoproduction in 
Ref. full . There the 77-photoproduction off the neutron, which 
is beyond the scope of this work, will also be investigated. 
The quality of our fits is given in terms of the^ of , 
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crN -6 V n(E) ^ 
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where z = cos(#), n(E) is the number of data points at energy E, 
cr is the number of distinct energies and N is the total number 
of data points. Moreover 5 = 4 denotes the number of degrees 
of freedom, which are called collectively ji. The letter e denotes 
the experimental values. Using now a random walk minimiza- 
tion procedure, we obtain a x\ f = 0.9997 with the following 
parameters 

logOWGeV) = -0.611, log(/v/GeV) = -0.512^, 
logOWGeV) = +1.845, logOWGeV) = -5.112^, 

where the error bars are obtained by varying the parameters 
such that the^ of is increased by one. No uncertainty is given, 
if it does not affect the value of a parameter at the given accu- 
racy. Some representative differential cross sections are shown 



in Fig. [2] Quite a good agreement between the model and ex- 
periment is achieved for energies more than 100 MeV above 
threshold, which can be seen from the plot of the integrated 
cross section in Fig. [3] Above the /^A threshold, the total cross 
section exceeds the data, this will eventually be overcome in a 
more precise NLO calculation. 

It is commonly believed that the first nucleon resonance, the 
S 11 (1535) saturates the cross section close to the tjN threshold. 
Moreover, it is known that this state can be understood as dy- 
namically generated already from the leading chiral order ver- 
tex, i.e. the WT-interaction. Since the latter is also the driving 
term of our hadronic amplitude it is worth to have a look at the 
modulus of the electric dipole amplitude Eo+ as a function of 
s = El ms on the second Riemann sheet. The resonance appears 
at 

(4) 



£ cms = (1525.9^-4 _ illl.4!j;jj MeV 



which is in good agreement with the extraction from phe- 
nomenological or other coupled-channel approaches collected 
in 111911 (at least for the real part) and our recent determina- 
tion from scattering data 1H2I1 . It can therefore be identified 
with the (1535) resonance, which is dynamically generated 
in the present approach. The imaginary part of the pole position 
is larger than the values of more recent phenomenological ap- 
proaches collected in |[l9ll . however the uncertainty given above 
reflects only the influence of the errors on the parameters and is 



3 





P?7 — tttTtt- 

/ 


AK 






***** 



1.50 1.55 1.60 1.65 

E cms [GeV] 



Figure 3: Total cross section for the best fit of our model (solid lines with shaded 
error bands) compared to the experimental data from the Crystal Ball at MAMI 
(filled circles with error bars). The vertical dashed lines denote the pr] and the 
KA thresholds. 

certainly underestimated. Moreover this large width is consis- 
tent with our determination from scattering at NLO. 
For the s-wave multipole Eq+ we obtain at threshold 

E 0+ = (-12.39!l;S + * 16 - 15 -L8s) ' 10" 3 /M, + . (5) 

Its modulus |Zso+| = 20. 35*23 J • 10~ 3 /M n + is somewhat larger 
than the one obtained in the early calculation in unitarized chiral 
perturbation theory [7] and in certain resonance models includ- 
ing the S ii (1535), see e.g. Ref. 12011 . Moreover, the ratio of the 
imaginary to the real part is about 1.06 .. . 1.68 and agrees with 
estimations from resonance models lEoll . 
Extension of the model: As stated before, one should include 
all NLO contributions in the hadronic as well as the electromag- 
netic part of the production amplitude. However, the full inclu- 
sion of such terms is beyond the scope of this work and will be 
given later As a first step to study the impact of such a 

modification, it is interesting to consider the magnetic moment 
couplings of the photon to the nucleon, since these play an im- 
portant role in phenomenological studies. The corresponding 
terms of the NLO chiral Lagrangian read 112111 

b 12 (Bo* y {F; v , B}) + b 13 (Bo* y [F; v , B]) , 

where F+ v = uF^u^ + u^F^u and F^ v = d^v v - d v v fi is the 
photon field strength tensor. The corresponding additional ver- 
tices are restricted to yB —> B and yBcp —> Bcp only, thus 
they do not give rise to any change of the hadronic part of our 
photoproduction amplitude. As a matter of fact unitarity in the 
subspace of meson-baryon channels is still preserved and the 
gauge invariance is guaranteed automatically by the functional 
form of the interaction itself. 

Within our model the inclusion of the £12,13 terms leads to 
an additional complication due to the A <-» iP transition in- 
duced by this terms. Without going into details, see [J^lllJ] for 
discussion, we have to modify the regularization scheme, such 
that the pure baryon integrals are set to zero from the begin- 
ning. Thus, our treatment of the loop integrals is, in effect, sim- 
ilar to the EOMS regularization scheme advocated in ref. i22ll . 
Moreover, the same transition implies that jjk = /^kh = Vka, 



which reduces the original parameter space of our model to 

We have performed fits to the Crystal Ball data on differen- 
tial cross sections in the same manner as before. However the 
energy region in question is reduced since no good agreement 
was achieved even for moderate energies. For the energy region 
1487. 8MeV < £ cms < 1541. 8MeV we can obtain a fit which 
minimizes the^ of to 1.75. For the same regularization scheme, 
but without the inclusion of £12,13 terms we obtain x\ of = 3.44. 
However, the agreement with the data is worse compared to the 
the leading order approach. The conclusion to be drawn is that 
the inclusion of the Pauli-terms ~ £12,13 alone does not improve 
the LO description. In fact, it is mandatory to perform a full 
NLO analysis. In particular, as the Pauli coupling is much for 
important for the neutron, we refrain from giving the LO results 
for yn — > rjn here. 

Summary and outlook: We have shown that the precise 
data of 77-photoproduction in the threshold region can be de- 
scribed rather accurately within the gauge-invariant chiral uni- 
tary framework. The S n(1535) resonance is generated dynam- 
ically and its pole position in the complex plane agrees with 
earlier determinations. To go to higher energies, to improve 
the precision and to investigate also the interesting properties 
of the same reaction on the neutron, one has to go to NLO and, 
in the long run, also include the Born terms in the unitarization 
scheme. Such efforts are under way Il4ll . 
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